Asymptotically exact solution of a local copper-oxide model 
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We present an asymptotically exact solution of a local copper-oxide model ab- 
stracted from the multi-band models. The phase diagram is obtained through the 
renormalization-group analysis of the partition function. In the strong coupling 
regime, we find an exactly solved line, which crosses the quantum critical point 
of the mixed valence regime separating two different Fermi-liquid (FL) phases. At 
this critical point, a many-particle resonance is formed near the chemical potential, 
and a marginal-FL spectrum can be derived for the spin and charge susceptibilities. 
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One of the most challenging theoretical issues has been raised by the experimental ob- 
servations on copper-oxide based metals: What is the appropriate description of the normal 
state which does not fit into the Fermi-liquid (FL) phenomenology. In an attempt to find 
the unifying features in the diverse observed anomalies, a phenomenological marginal-FL 
(MFL) spectrum for the spin and charge fluctuations was proposed by Varma et al. [1]. 
The essential point is that the frequency dependence of the susceptibilities is singular, while 
the momentum dependence is assumed smooth. It is expected that the study of the related 
impurity models would shed some light on the basic physics of the lattice models, in the 
same sense as the Anderson model vs the single band Hubbard model. 

Recently, Varma et al. [2,3] have explored the local non-FL properties of a generahzed 
Anderson model, equivalent to a single-impurity version of one of the multi-band models 
proposed to describe the physics of the copper-oxide compounds [4,5]. In fact, such a model 
was first studied by Haldane in the late 70s on Anderson's suggestion [6]. They reahzed 
that the usual Anderson model should be supplemented by screening channels to saturate 
the Friedel sum rule, which is crucial for the understanding of the mixed valence physics. 
The recent Wilson's renormahzation-group (RG) study [2] and analytical arguments [3] 
have provided a microscopic scenario for the MFL phenomenology within this generalized 
Anderson model. 

In this Letter, we derive explicitly physical properties of this local copper-oxide model 
through exactly solving the Hamiltonian in the asymptotic limit. We first separate the 
model into "hybridizing" and "screening" parts. The first part is the usual Anderson model 
plus X-ray edge-like (XRE) scattering terms, which were neglected in all previous studies 
[2,3,6]. This part reduces to a hybridization process and potential scatterings, which can 
be merged into an effective hybridization via a canonical transformation, giving rise to a 
local FL behavior. The second part is the multi-channel XRE scattering terms, and can 
be reduced to a single spinless channel, exhibiting the Anderson catastrophe. These XRE 
singularities can also be transformed into the effective hybridization, so the physics of the 
full model results from a competition between the above two different physical factors in the 
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hybridization. Then we apply the RG theory [7] to derive flow equations from which a phase 
diagram is found. There exist two different FL phases: a Kondo phase and a local free FL 
phase separated by a mixed valence phase, displaying a local MFL behavior. Finally, from 
the partition function and the effective Hamiltonian, we take the strong coupling limit, and 
an exactly solved hne (analogue of the Toulouse hmit in the Kondo problem) can be found, 
which crosses the quantum critical point of the mixed valence phase where a many-particle 
resonance is formed between the localized impurity and the conduction electrons. Properties 
of the above three phases, especially, the mixed valence phase controlled by the critical point 
can be exactly derived. 

The local copper-oxide model is specified by the Hamiltonian [2] 

H — Hh + Hg, 

Hh = Y. ^kCla,oCk,a,o + erf^d + Und,^nd,i + J^i^l^fida + h.c.) 

+ H -^Ck,a,oCk',afl{nd,a - ^ + X! -^Ck,<T,oCk' ,<j,o{nd,a - 1^ , 

Hs^ (^kCl^a,lCk,a,l + -^^k,a,fik',<j,l{'^d- -^)-, (l) 

k,(j,l>Q k,k',a,l>0 ^ 

where Ud — Yl,a '>T'd,a, N is the number of the lattice sites, and we have separated the Hamilto- 
nian into hybridizing and screening parts and distinguished parallel-spin and opposite-spin 
XRE scatterings in the hybridizing channel. The localized impurity hybridizes only with 
channel I — 0. The chemical potential /x is set to zero, and we are interested in the case 
when the local impurity level is close to zero. The spinless version of this model, i.e. the 
multi-channel resonant level model, has been solved exactly [8,9], displaying a FL vs non-FL 
transition as the interaction parameters are varied. In fact, is the usual Anderson model 
plus XRE potential scatterings, while Hs is the multi-channel XRE Hamiltonian. It is more 
convenient to take an infinite U limit as in the usual treatments for the Anderson model (the 
low-energy physics is kept). Thus, we add a local constraint for the local impurity < 1. 

First, we use abelian bosonization to handle the XRE singularities of the screening 
channels, which reduce to a one-dimensional problem with only one Fermi point for each 
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channel and the dispersion is hnearized with a cutoff ko [10]. Since the spin degrees of 
freedom of the screening channel electrons are trivially involved, they can be separated from 
Hg. Moreover, we can assume Vi = Vg for alH > without loss of generality, so the channel 
index can be dropped. Thus, all the screening channels are described by a single spinless 
channel [3]. The resulting form is 

H's = E ^4«. + E ^s\l^{ai + a,){n, - h, (2) 

A;>0 P k>0 * ^ 

with Vs = y/2NsVs and Ng is the number of the screening channels. a\ and are the bosonic 
operators describing the charge degrees of freedom of the screening channel and p = (hvp)'^ 
is the density of states at the Fermi point. Employing the inverse bosonization, we can 
transform bosons back to fermions: = J2k ^kslsk + jf J2k,k' sl.Sk'{nd—^). The Hamiltonian 
(0) can be diagonalized through a canonical transformation [8]: U = exp{Y,k>o ^/^{(^k ~ 
a\){nd — |)}, where 5s ~ T^pVs is the phase shift generated by the XRE scattering of all 
screening channels at the impurity. 

For the hybridizing channel, employing another canonical transformation [8,11]: S = 

exp{E{k>o,a} ^{h,a - bl^^){nd,a - ^)}exp{J2{k>o,a} ^{bk,a-bl„){nd,a-l)}, where 6o ^ 
tt/jVo and 6q ~ TrpFj are the phase shifts induced by the parallel-spin (cr) and opposite-spin 
(a) XRE scattering of the hybridizing channel, respectively, we merge XRE scatterings into 
the hybridization and obtain 

k k 
Heff= E -b\,abk,a + ^'d'nd + tY^{/^i)ld„ + h.c.) + Y,-(^W, (3) 

fc>0,(T P <y k>0 P 

with i/j^ = Vk^exp[j:k>o^^^ibk,a - bl^)]exp[-Ek>o ^ibk,a - &U)], A = 
exp[Y,k>o ^/^i'^k — a|.)], e'd = — pVqVq, and < 1. A fermionic form is obtained by 
inverse bosonization 

Heff = E ^kClAa + e'dUd + tJ2i^^ld^ + h.c.) + E ^ksW, (4) 

k,cr o" k 

where = ^ = -^TkCk,., A = (M~^(^o)^, ^^o = 

^EfcSfc, and Ud < 1. 
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Next, we derive the partition function, dividing Hejf as Heff = Hq + Hi with Hq and 
Hj as the free and the hybridization parts of (^, respectively. Parallehng all the strategies 
of previous studies [6-8], we write the partition function in terms of a sum over histories 
of the impurity. Each history is a sequence of transitions between the three local d states 
I a >=\ >, and \ a >, a =t,i- The transitions take place at the imaginary time 
0<ri<...<r„</?=|;: along the Feymann trajectory the local state is at | (Tj+i > from 
Tj to Tj+i {i=l to n). The partition function of (||) is now given by 



-0 {o-,(Ti,..,(T„+i=cr} {a=0,cr} 



TT TT TT T 

where the bare hybridization strengths are defined as F = pt^, while the cutoff factor r = 
The effective "magnetic field" reflects the differences of the local state energies: i?^ = — e'^, 
-^d = e'^. The fugacity (2/". ^.^^ vTV) is the amplitude associated with a transition from 

I ai > to I a^+i >, with = ^.^...^.-i, ?/°2,+i,<x2. = 1; Vl^.^o^.^^ = h V^.^+^a,, = 

5o-2i+i,(T2,- The long-range logarithmic interaction between the flipping events arises from the 
reaction of the conduction electron bath towards the transition between the local states. The 
local disturbance on the bath involves two factors: the absorption or emission of the local 
conduction electrons and the change in the local potential that the conduction electrons 
experience [8]. Both kinds of disturbance are incorporated in the effective "charge" factor, 
i.e. the coefficient of the logarithmic function of (^. 

The partition function could be obtained directly from the model Hamiltonian (|1]) 
without bosonization, using the famous fermion techniques [12] in certain asymptotic limit 
which, we believe, is also valid here. This alternative derivation would allow us to rec- 
tify the phase shifts obtained by the bosonization treatments to the exact expressions: 
5s = 2tan^^{^pVs), So = 2tan^^{^pVo), and S'q = 2tan^^{^pV'o) so that the corresponding 
effective Hamiltonian @ might be used beyond the range of validity for the bosonization 
method, especially in the following strong coupling limit where the renormalized parameters 
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of the model recover their bare values. 

To set up the RG flow equations, we can directly employ the scaling theory proposed by 
Anderson et al. [7] in the Coulomb gas representation. The RG equations describe the flow 
behavior as the bandwidth is reduced. They are given by 

dr ^ dcd r (^7 , 

= -7r, -jT^^-, 37rz«^-(7 + irrr, (6) 



dlnr dlnr n dlnr 

with 7 = -^ + (^)^ + (^)^ + (^)^ describing the total interaction strength between the 
conduction electrons and the local impurity, which should be positive in our case. These 
equations were derived by assuming Fr < 1, a rare gas of spin-flips. In the zeroth order, we 
can construct two invariants (r*T* = 1) : T* — r(rro)i^ and — — (^'Tq)^], 

where €d, F, and Tq are initial (bare) parameters. In terms of these scaling invariants, 
the running resonance width and impurity level are written as T{t) = (r*y~"' [t)~'^ and 
ed(T) = + ^[1 - (r*r)-T]. Obviously, for the case = = K = 0, i.e. a trivial 
hmit of 7 = 0, the two expressions become r(T) — F* and ed{r) — e*^ + ^/n(r*T), which 
exactly recovers the Haldane's RG results for the standard Anderson model [7]. Moreover, 
a complete phase diagram can be determined by comparing the invariants and F*. In the 
plane F* — 7 (Fig.l), there are three phases corresponding to different impurity-occupancies: 
single-occupancy regime (e^ <C — F*) where a singlet state is formed and the Kondo effect 
shows up; zero-occupancy regime (e^ » F*) where the model corresponds to a local free 
FL phase, and the mixed valence regime (| |< F*) where < > fluctuates between 
and 1 phases. The mixed valence regime separates the < Ud >= 1 and the < >= 
phases with crossover hues F* Ri and F* ^ j^::^: respectively. Although parameter 7 
is not renormalized in the zeroth order of Ft, in its flrst order it is renormalized to smaller 
values as r increases. On the other hand, from F* pa —jTre^ + F, F* increases when > 
and decreases, while it decreases if < and decreases in absolute value. Hence the flow 
directions indicated in Fig.l. In the end, the renormalized 7 tends to zero as — > 0, F* — > F. 
However, when F(r)r ~ 1, the RG equations are not correct quantitatively. 

In addition to this RG analysis, a strong coupling limit can be extracted independently 
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from the partition function (5) and the effective Hamiltonian (4). As follows from (4), when 
opposite-spin XRE scattering in the hybridizing channel renormalizes to zero while the XRE 
scatterings of the parallel-spin and the screening channel reach their respective unitary limit, 
i.e. 6o = Ss = and (5q = or 7 = 0, the phase shifts due to hybridization and parallel- 
spin XRE scattering in the hybridizing channel compensate each other, thus the hybridizing 
electrons become completely free. In such a strong coupling limit, the effective Hamiltonian 
becomes 

Ht = J2 ^kClaCk,a + £dnd + tJ2isUi + h.c.) + J2 ^ks\sk, (7) 

k,cr <y k 

with constraint rid< \ oy n^.^ SqSq = 1, which reflects the Friedel sum rule in this limit. It 
is obvious that the partition function derived from the Hamiltonian (|^) is exactly the same 
as (H) for = 8s = and 5g = 0, z.e. 7 = 0. In this sense, the strong coupling limit found 
here is somehow analogous to the Toulouse limit of the Kondo problem [13], although the 
actual physics involved is quite different. The most essential difference is that the unitary 
limit has been actually reached in our case. The vanishing of opposite-spin XRE scattering 
in the hybridizing channel is exactly what is required by the infinite U limit, because any 
finite hybridization between opposite-spin hybridizing electrons and the local impurity (to 
compensate XRE scattering) will contradict the single occupancy constraint. 

In Fig.l, 7 = is a strong coupling limit line of this local copper-oxide model. In (^, 
only charge of the local impurity a = + d\) is coupled to the conduction electrons, 

while the spin /3 = — is decoupled except for the constraint. Thus, in this limit, the 
Hamiltonian is: Ht = Hk ^k^lsk + edii^-a + ^/3) + V^t{sla^ + h.c), where Ua + np <1 and the 
hybridizing electrons do not show up explicitly. This Hamiltonian is essentially the same as 
Eq. (14) of [3]. Since Ht conserves rip, we can calculate physical quantities by taking the 
trace on its two subspaces = 0, 1. 

(i). When ^ edc, a critical value to be defined later, the tIq, = = state is favored 
in the low-energy regime. All charge fluctuation processes are frozen out and SqSq = 1 at 
the impurity site so that the Friedel sum rule is saturated. The hybridization strength F* 
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should be renormalized to zero, and a local free FL behavior is thus displayed [14]. 

(ii) . The opposite case <^ e^c favors Ua+rip = 1 in the low-energy regime, and there are 
two possible configurations: = 0, n/3 = 1; Ua = 1, n^j = 0. All charge fluctuation processes 
are also frozen out but SqSq — 0. The hybridization strength F* should be renormalized to 
+00, and the system is scaled to Wilson's strong couphng fixed point of the Kondo problem: 
a local FL behavior in its unitarity limit [14]. The symmetry of the ground state for the 
present case is different from that at ^ e^c, so we anticipate a quantum critical point at 

(iii) . When edc, the localized impurity fluctuates between zero- and single- 
occupancy < ria + >— > |, and hybridizes with only part of the screening electron 

< SqSq >— > |, corresponding to the mixed valence phase. At the special point (7 = 0, 
F* = F) along the strong coupling line, the above two different FL states are degenerate. 
This special point is just the quantum critical point controlling the physics of the whole 
mixed valence phase. An analogous quantum critical point was found in the two-channel 
or two-impurity Kondo problems [15]. At zero temperature, we find e^c ~ _31^y and 

< Ha >~ \, < Up >Ri 0. Thus, the local impurity level is close to the chemical potential. 
Using the phase-shift representation of the Friedel sum rule: < > +^Sh{lj) + :^(^s(a*) = 1, 
we easily obtain the phase shift of the screening electrons caused by the final hybridization 
is I at the chemical potential. Since both hybridizing and screening electrons are involved 
and there is a constraint on the local impurity, the hybridization becomes a many-particle 
resonance, drawing some weight of the one-particle spectra from higher energies at the scale 
of the charge transfer gap in the insulating state. Such a many-particle resonance breaks 
down the Landau correspondence between the low-lying excitations of the interacting and 
non-interacting fermions [3]. At finite low temperatures, the impurity charge and longitu- 
dinal spin susceptibilities Xa,p can also be calculated using the relations — {a^P + /S^a), 
p — {a^a + P^P ~ I)- If hs-s been found that Xa,p are proportional to F~^/n(^), as expected 
from the MFL phenomenology [3] . Thus, the MFL behavior controls the whole mixed valence 
regime. 
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In conclusion, we have presented an asymptotically exact solution of the local copper- 
oxide model including the charge fluctuations of the screening electrons. The physical picture 
of the breakdown of the FL behavior pointed out in Ref.[3] is basically correct. However, 
the justiflcation of their physical arguments involves several unclear approximations. The 
crucial point is that in all the previous studies [2,3,6], the parallel-spin XRE scattering in the 
hybridizing channel was assumed to be zero, i.e. Vq = 0. However, namely this assumption 
obscures the physical features of the model and makes the problem much more involved. 
Due to the emergence of new relevant variables in their theory, in principle, they can not 
reach the strong-coupling limit. Since our solution is based on both RG flow analysis and 
the strong coupling effective Hamiltonian, the MFL behavior should be a universal property 
of the mixed valence phase. Of course, whether a speciflc system is at the quantum critical 
point depends on a special combination of parameters. A more interesting question, whether 
the chemical potential of a real mixed valence system is pinned at the critical point, requries 
further studies. 
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Figure Caption 

Fig.l. The phase diagram of the local copper-oxide model. I, II, and III are single- 
occupancy phase, mixed valence phase, and zero-occupancy phase, respectively. The two 
thick lines correspond to the crossover lines, while the two slender lines with arrows roughly 
show up the flow directions and the horizontal one is = line. 
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